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The post-Newtonian expansion appears to be a relevant tool for predicting the gravitational 
waveforms generated by some astrophysical systems such as binaries. In particular, inspiralling 
compact binaries are well-modelled by a system of two point-particles moving on a quasi-circular 
orbit whose decay by emission of gravitational radiation is described by a post-Newtonian expansion. 
In this paper we summarize the basics of the computation by means of a series of multipole moments 
of the exterior field generated by an isolated source in the post-Newtonian approximation. This 
computation relies on an ansatz of matching the exterior multipolar field to the inner field of a 
slowly-moving source. The formalism can be applied to point-particles at the price of a further 
ansatz, that the infinite self-field of point-particles can be regularized in a certain way. As it turns 
' out, the concept of point-particle requires a precise definition in high post-Newtonian approximations 

of general relativity. 
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I. INTRODUCTION 



Binary systems of compact objects (neutron stars or black holes) emitting gravitational waves in their late stage of 
I \ evolution leading to a final coalescence might play the major role in the XXIth-century gravitational-wave astronomy. 
^ ■ In the final stage, the two objects collide; for instance two black holes form a single black hole which rings down and 
C*~) ' emits gravitational waves before settling down into a stationary configuration. We refer to the papers in this volume 
devoted to numerical relativity |Q and to the close-limit approximation || for descriptions of the coalescence phase. 

In the earlier phase, preceding immediately the coalescence, the two objects undergo a long (adiabatic) inspiral 
driven by the emission of gravitational radiation, or equivalently by the radiation reaction forces applied to the orbit. 
In this paper we are mainly interested in the inspiral phase of compact binaries. During this phase the gravitational 
radiation is essentially produced by the dynamical motion of the two compact objects. In principle, the dynamics 
can be well-approximated by a post-Newtonian expansion of general relativity. Given the highly relativistic nature 
of inspiralling compact binaries (the orbital velocity can reach 30% of the speed of light in the last rotations), the 
problem is just that of pushing the post-Newtonian approximation farther enough in order to be useful to future 
observations. In recent years it has been realized by several groups |3]-|To[| that in the case of inspiralling neutron 
star binaries the post-Newtonian expansion should be controlled up to the very high 3PN order. During the inspiral, 
the internal structure of the stars plays a little role, and one can conveniently describe the two compact bodies by 
"point-particles" (what only matters are the two masses, but apart from that the objects can be ordinary neutron 
; I ■ stars, or black holes or even naked singularities). Note that a priori the concept of point-particle does not make sense 
in general relativity except in the test-mass limit. However within a post-Newtonian approximation one can give a 
sense to what we call a point-particle; but even there, particularly when going to high post-Newtonian orders, the 
concept of a point-particle is non-trivial and must be carefully defined. 

A post-Newtonian computation of the inspiral of compact binaries can be based on the following strategy. One 
starts by implementing a general formalism for the dynamics and the gravitational-wave emission of a slowly-moving 
isolated source. By slowly-moving we mean the existence of a small post-Newtonian parameter, say 
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where T^ v denotes the matter stress-energy tensor in some Cartesian coordinate system covering the source; Q Greek 
indices take values 0,1,2,3, and Latin 1,2,3. We allow e to be as large as 30% in order to cover the case of inspiralling 



*In Black Holes and Gravitational Waves, Proceedings of the 9th Yukawa International Seminar, T. Nakamura and H. Kodama 
(eds.), Prog. Theor. Phys. Suppl. No. 136, 146 (1999). 
1 We assume that the source is self-gravitating, so that e ~ (GM/ac 2 ) 1 ^ 2 where M and a denote the mass and radius. 



1 



compact binaries. As we see, strictly speaking the source is "slowly-moving", but in fact it can be quite relativistic. As 
usual, neglecting terms in the gravitational field according to their formal order in e is valid in the so-called near-zone 
of the source; the region for which r/X = 0(e), where A is a typical wavelength of the radiation field. 

The crucial demand on any wave-generation formalism is to be able to relate the gravitational waveform far from 
the source (to order l/R in the distance to the source at retarded times) to the matter stress-energy tensor T^ v . This 
problem of relating the retarded far-field to the source's matter content is extremely difficult within the exact theory, 
because of the non-linearities of the field equations; however, the solution exists in a framework of post-Newtonian 
approximations for slowly- moving sources. Another requirement for a general formalism is for its ability to control the 
equations of motion of the source, and in particular the gravitational radiation reaction forces therein. The radiation 
reaction forces are to be consistent with the radiation field at infinity (by definition) , and thus to depend on boundary 
conditions such as the no-incoming radiation condition imposed at past null infinity (which ensures that the source 
is physically isolated). In a post-Newtonian approximation, a difficulty of the problem is that the radiation reaction 
forces enter the equations of motion of the source which are determined by the post-Newtonian expansion valid only 
in the near zone. Therefore one must supplement the post-Newtonian expansion by a condition of matching the near- 
zone field to the radiation field [see (2-5) below]. In this article we outline a particular post-Newtonian formalism, 
making extensive use of multipole moments, which is issued from work of Blanchet and Damour, JH] Damour and 
Iyer, JL4| and Blanchet. fi6H18|l A different formalism recently defined by Will and Wiseman |H| on foundations laid 
by Epstein and Wagoner J20[ and Thorne is described by C. Will in this volume. 

Note that the present post-Newtonian formalism is interested in the formal post-Newtonian expansion e — * 0; in 
particular it does not try to investigate the exact mathematical nature of the post-Newtonian series. Thus the con- 
nection between the approximation and the exact theory is not controlled^] Simply, one looks for a mathematically 
well-defined scheme for generating the successive post-Newtonian approximations. The hope is that in applications, 
when comparing to observations the post-Newtonian prediction developed to high enough order in e, we shall drasti- 
cally limit the magnitude of systematic errors, resulting in a very accurate measurement. 

Be aware that at present no systematic algorithm for generating the post-Newtonian series of the near-zone gravita- 
tional field exists. One can express the multipole moments of the source in terms of the infinite formal post-Newtonian 
series, Jl8[ but one does not know how to compute explicitly this series to all orders, in the manner, say, of the ex- 
plicit algorithm for generating the post- Minkowskian series in the source's exterior. |flfl B ut of course, the first few 
post-Newtonian approximations, up to, say, the 2.5PN order, are very well understood. [|23|-p8| A difficulty associated 
with the post-Newtonian iteration is the appearance in higher post- Newtonian orders of Poisson-type integrals which 
are divergent at the bound at infinity. This comes from the fact that the post-Newtonian iteration is valid only in the 
near-zone, and so the post-Newtonian coefficients rapidly blow up far from the system. This signals that the standard 
Poisson integral does not constitute the relevant solution of the Poisson equation in this context. 

A crucial element of the present formalism is that it is a priori only valid for continuous (in fact, smooth) matter 
distributions. Thus, T^ v is assumed from the start to be regular, for instance to describe a smooth hydrodynamical 
fluid. This excludes a priori the very interesting application to "self-gravitating" point-particles (as opposed to "test" 
point-particles), whose self-field becomes infinite at the location of the particle, thus creating a singularity So we 
are obliged to introduce a new ingredient. Our proposal is that this be a regularization a la Hadamard p9| , [30| for 
removing systematically the infinite self-field of the particles. Only when we assume a regularization, can we use 
a T^ v constituted of Dirac delta- functions. The removing of the divergent terms is for the moment done without 
further justification; simply this is an ansatz which, as far as we can see, yields consistent computations in practice, 
and which has been checked to yield the correct result in some cases, but which remains an ansatz. Notably we do not 
prove that the regularization is still permissible in higher post-Newtonian approximations, or even that it is possible 
to find a consistent regularization at all, or that the result of two different (though consistent) regularizations would 
be the same. We will content ourselves with the use of the Hadamard regularization which yields in practice, to the 
first few post-Newtonian orders, some consistent (and rather elegant) computations. 

A different strategy is possible when we have at our disposal a natural background space-time. This is the case 
when the mass ratio of two particles is so small that one can view one of them as moving in the Schwarzschild or Kerr 
background generated by the other. In the small mass-ratio limit the particle moves on a geodesic of the background, 
and we can compute the emitted radiation using a linear background perturbation. This approach has reached a 
mature state: notably the gravitational radiation emitted by a test particle in orbit around a Schwarzschild black 
hole was computed to very high post-Newtonian order; [pl|-|33f this computation represents an important benchmark 
against which the standard post-Newtonian expansion can be checked. To second-order in the black-hole perturbation 



2 A mean to understand the limit relation of Einstein's theory to Newton's is to introduce a frame theory a la Ehlers. [G2 
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there are terms vanishing in the mass-ratio limit, the particle reacts to the gravitational radiation emission, and again 
we face the problem of how to regularize the particle's self-field. 



II. GENERAL SOURCES 

For a general compact-support stress-energy tensor T^ v , we want to solve the field equations of general relativity in 
the form of a post-Newtonian expansion e — > 0. We reduce the field equations by means of the condition of harmonic 
coordinates, that is d v h^ v — where W v = yf—gg 11 " — r\^ v , so that 

Uh^ = —^\ g \T^ v + A^(h, dh, d 2 h) , (2) 

where we have introduced a Minkowskian background, rf v — diag(— 1, 1, 1, 1) = 1}^ (with respect to which all indices 
are raised and lowered), and the associated d'Alembertian operator, □ = □,, = rf v d ii d v . The gravitational source 
term A'*" is a complicated functional of the field and its first and second space-time derivatives. We define the total 
stress-energy (pseudo-) tensor r Miy of the matter and gravitational fields as 

T ,u = i | TM „ + J?—K^ . (3) 
lonG 

Of course t^ v is not a generally-covariant tensor, but it is a Lorentz tensor relatively to our Minkowskian background. 
It is conserved in a Lorentz-covariant sense, and this is equivalent to the covariant conservation of T^ u , 

d v T^ = V„T^ = . (4) 

The propagation of h^ v subject to the Einstein field equations (2-1) is a well-posed problem, for which we need to 
choose some initial conditions in the past. In this paper we shall assume that the field is stationary before some 
remote date — T, so that there is no radiation generated by sources at infinity incoming onto the system. Arguably, 
the condition of stationarity in the past is too strong; for instance it does not cover a physical situation where two 
bodies moving initially on unbound (hyperbolic-like) orbits would form a bound system. However one may in some 
cases justify the assumption a posteriori, by checking that the formulas obtained under it are still valid in a more 
general physical situation such as the initial scattering of two bodies. With no-incoming radiation one can transform 
the differential Einstein equations (2-1) into the integro-differential equations 

r= « D - V ^-^ / J^ T ^( X ',t-\ X - X '\/c) , (5) 

where D^ 1 denotes the standard retarded inverse dAlembertian. 

In our approach we resolve the wave-generation problem by finding the solution of the field equations (2-4) in the 
exterior region of the source (outside the compact support of T^) by means of an infinite multipolc-moment series 
for W, that we denote M.{h^ u ). This means in particular that we must be able to relate the multipole moments 
parametrizing this series to the matter content of the source. In general this is not an easy task, but this can be 
done using the post-Newtonian expansion, in the physical case where the source is slowly-moving. The latter source 
multipole moments are then "propagated" with the help of a post-Minkowskian expansion to large distances from the 
source, and related there to the so-called radiative multipole moments directly accessible to a far-away observer.^] 

To obtain the expression of the multipole moments of the source in terms of T^ v we use an asymptotic matching 
between the multipole expansion M-Qi 1 ^"), valid everywhere outside the source, and the post-Newtonian expansion 
denoted h , valid in the near-zone. For slowly-moving sources the two domains of validity of the multipole and 
post-Newtonian expansions overlap in the so-called exterior near-zone, and we can write there the numerical equality 
M{U LV ) = \i . Then we transform this equality into a "matching equation" , that is an equation between two series 
of the same nature, which is formally valid "everywhere" . For this purpose we replace the multipole terms on the 
left-hand side of the equality by their formal post-Newtonian expansions (this means expanding all retardations t—r/c 
when c — > oo; so in fact the latter expansion is equivalent to a near-zone expansion r — > 0), and we consider on the 



3 Using a post-Minkowskian expansion for the exterior field in conjunction with the multipolar series is an old idea of Bonnor, 
B4| later generalized by Thorne. pif 
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right-hand side the multipole decomposition of (each of the coefficients of) the near-zone expansion. Therefore the 
matching equation reads 



M(h^) = M(hT) . (6) 

Satisfying this equation permits in principle to find the unique physical solution of the field equations valid inside and 
outside the source (and, in the present case, in the uniquely defined harmonic coordinate system). 

Let us simply announce here the result |k|[18| for the expression of M(h^ v ) satisfying at once the vacuum Einstein 
field equations (with no-incoming radiation) and the matching equation (2-5). There are two terms, 

M(hn = --rE hr dL r^ {t r,c) + FPs=0 d « ^(A"")] . (7) 

The second term constitutes a particular solution of the vacuum field equations, which is defined by analytic continua- 
tion in a complex parameter B as being the finite part of the Laurent expansion when B — > (in our notation FPs = o). 
The reason for the need of a finite part is that the integrand of the retarded integral is in the form of a multipole 
expansion which is not valid inside the source, and which is actually singular at the spatial origin of the coordinates 
r = located inside the source. Now the first term in (2-6) (in which we denote L = i\ ■ ■ ■ i; and 8l = 9^ ■ ■ ■ d^) 
represents clearly an homogeneous solution of the wave equation. This solution is uniquely singled out so that when 
added to the particular solution it ensures the satisfaction of the matching equation (2-5). The "multipole moments" 
parametrizing this solution are given as |l8| 

T^{u)=¥V B=0 J d 3 x \x\ B x L J dzStWr^fau + zlxl/c), (8) 

(where u = t — r/c). The moments are also defined by analytic continuation in B, which makes them well-defined 
mathematically (note that the integrand behaves typically as a positive power of the distance at infinity, and thus 
that the integral would be strongly divergent at infinity without any finite part). In (2-7) we use a special notation 
for a symmetric trace-free product of vectors: ij, = STF(xl) where xl = Xi t • • -x^ (we have 6i l i ll XL = 0). The 
z-integration in (2-7) involves the weighting function 



*»(*) = 2 i+H\ {1 ~ Z > ' / <feW = l- (9) 

In the limit of large I we have lim^oo Si — S the Dirac measure. 

The crucial point about the multipole moments (2-7) is that they are generated by the post-Newtonian expansion of 
the stress-energy pseudo-tensor, that is t mi/ , rather than t^ v itself. It is at this point that our assumption of matching 
with a slowly-moving post- Newtonian source enters. Note that although the integrand of the multipole moments is in 
the form of a post-Newtonian expansion valid only in the near-zone, the integration is to be performed on the whole 
3-dimensional space (see Ref. Eq| for details). This is one of the beauties of the analytic continuation, that it permits 
to handle integrals over M 3 without introducing a cutoff at the edge of the near-zone. On the contrary, Will and 
Wiseman [ fl9| do not use analytic continuation, and consider integrals extending only over the near-zone. It can be 
shown that the present formalism is equivalent to that of Will and Wiseman. 

With the expression (2-7) in hands, it is straightforward to define six sets of irreducible multipole moments 
{II, Jl, Wl, Xl, Yl, associated with the six independent components of T 1 ^ (ten minus four because of the har- 
monic coordinate condition). This constitutes our primary definition for the source multipole moments. With more 
work (considering a coordinate transformation in the exterior zone) we can further define only two sets of irreducible 
moments {Ml, Sl,} which are less directly connected to the source but are very useful in practical computations. 

Convenient notions of the source multipole moments being chosen, let us relate them to the radiative moments 
at infinity. We follow the standard definition |2f| l that the radiative moments parametrize in radiative (Bondi-type) 
coordinates the leading term 1/R in the distance to the source. Clearly, looking at (2-6), we see that the radiative 
moments have two contributions. Essentially the contribution coming from the first term of (2-6) is that of the source 
multipole moments II, J hi " • ' 1 ^h- In a linear theory the radiative moments would contain only this contribution, i.e., 
they would agree with the source moments. But in general relativity we have also the contribution of the second term 
in (2-6) (which is at least quadratic in the field strength G), so there will be many non-linear interactions between the 
source moments to any order in G. Therefore the radiative moments are given by some (very complicated) non-linear 



functionals of the source moments. We compute these with the help of a post-Minkowskian algorithm, llflq] That 
is, we are able to rewrite (2-6) in the form of a formal expansion when G — * 0, 
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M{hT) = Gh1 v [I L , J L , ■■■} + G 2 h^ + ■■■ + G n h^ + ■ • • , 



(10) 



where the first term is a solution of the linearized field equations which depends on the source moments [this is 
essentially the first term in (2-6)], and where all the subsequent non- linear corrections are constructed by post- 
Minkowskian iteration (see Ref. |ll]] for the proof that this can be done to any post-Minkowskian order). Next we 
change coordinates from harmonic to radiative ones (T, X) and expand the metric when R — > oo. All the physical 
information about the radiation field is contained into the so-called transverse-traceless (TT) projection of the spatial 
(ij) metric. The radiative moments {Ul,Vl} are defined from the term of order 1/R as 



X(^) TT (X,T) = 




where N{ = X l /R, Nl-2 — N n ■ ■ ■ Ni t _ 2 , etc., and where the TT projector reads Vij a b — (S ia — NiN a )(8 3b — N 3 N b ) — 
\{5ij — NiNj)(S a b — N a Nb). Now it remains to compare the result of the iteration of the field (2-9) with the relation 
(2-10) in order to deduce the radiative moments Ul and Vl in terms of the source moments, and therefore in terms 
of via the post-Newtonian expansion of the pseudo-tensor t^ v [see (2-7)]. 

The wave-generation is complete in this way. Note that the formalism is general in the sense that we did not specify 
a particular form of T^ v . Thus, for application to a particular problem, we need first to choose a model of T M1/ , for 
instance a perfect fluid or a model of point-particles (see below), and insert that model into the formalism. 



III. POINT-PARTICLES 



The choice of a model of point-particles in (post-Newtonian) general relativity is non-trivial because it is intimately 
related with the choice of a regularization for removing the self-field of point-particles. Of course, for a test particle 
moving on a fixed smooth gravitational background g^ v , there is no problem and the stress-energy tensor reads 

TZ^t) = - , (12) 

s j-g%h t )Vv"l<? V~9 B (x) 

where = dy^ / dt is the particle's coordinate velocity and S is the 3-dimensional Dirac measure (a particular case of 
distribution) . 

However, in the case of " self- gravitating" particles contributing to the gravitational field, the previous stress-energy 
tensor does not make sense. Indeed, already at the Newtonian order, the metric g pL , generated by two particles for 
instance contains the Newtonian potential U = Gm\jr\ + Gmijri where r± t 2 = |x — 2 | is the distance between 
the field point and the source points. Clearly U is infinite at the location of the particles, and must be regularized in 
some way in order to speak of its value at point 1, say {U)%. We know that the correct result in Newtonian gravity 
is (U)i = Gm,2/ri2 where r\i = \yi — y 2 \- But in high post-Newtonian orders we shall meet more difficult quantities 
such as (C/ 4 )i which necessitate a precise definition. 

Let us model the stress-energy tensor of two point-particles in post-Newtonian approximations of general relativity 
by 

T^ v (x, t) = — 1 1 v (13) 

V-(5p<t)i«/c 2 y/-g(x) 

where 1 <-» 2 means the exchange of 1 to 2 on the two previous terms and, {g pa )\ denotes the regularized value of the 
metric g pa at the location of the particle 1, according to a procedure based on the Hadamard [^9|,^0) partie finie of 
a singular function and a divergent integral. Note that in fact the chosen regularization defines our model of point- 
particles. A conjecture would be that another type of regularization, if sufficiently powerful to produce unambiguous 
results at high post-Newtonian orders (as seems to be the case of the Hadamard regularization), would yield identical 
physical results. 

To define the Hadamard regularization we introduce an appropriate class of functions T on 1R 3 . By definition F(x) 
belongs to T if and only if it is smooth on M 3 except at two singular points y 1 and y 2 , around which it admits a 
power-like expansion of the type 
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i? = E r ?/i(«)( n i) whenn^O (14) 

a 

(and idem for 2). The summation index a is assumed to be bounded from below, a > —ciq (no essential singularity), 
and to take real discrete values, say a € {aj}jgjv where dj £ JR. In fact in most cases one can assume a £ Zi. As 
indicated the coefficients /i( a ) depend on the unit direction ri\ of approach to the singularity. In (3-3) we do not write 
any remainder for the expansion because we shall need only formulas depending on the /i( ), and never the complete 
expansion itself. ^ 

First we define the Hadamard partie finie of the singular function F at the location of the singularity 1 for instance, 
say {F)±. Simply we pick up the coefficient of the zeroth power of r% in (3-3), namely /i(o), and average over all 
directions n\. 



(F)r ^ / ^/ 1(0) . (15) 

Second we define the partie finie (in short Pf) of the divergent integral J d 3 x F. We remove from M 3 two spherical 
balls surrounding the two singularities, of the form n,2 < s, where s is a small radius. Using (3-3) it is easy to 
determine the expansion when s — > of the integral extending on M 3 deprived from these two balls, next to subtract 
all the terms which are divergent when s = 0, and then to take the limit s — > of what remains. The result is the 
Hadamard partie finie. Actually we find some logarithms of s in the expansion, so there is an ambiguity linked with 
the freedom of choosing a unit of length to measure s. To account with this ambiguity we introduce a constant length 
scale si (and similarly S2). We find 

a+3 



Pf J d 3 x F = lim{/ d 3 xF+ J2 £3 / dUiAc 

J yJ t 2 >s a+3<-l J 



+ ln^J J rfr> 1 / 1( _ 3) + 1^2| , (16) 

The two definitions (3-4) and (3-5) are closely related to each other. To see this, apply (3-5) to the case where the 
function is actually a gradient diF (it is clear that F e T implies diF E T). We find 

Pf J d 3 x d l F = -4n(n\r%F)i - 47r(n 2 r|F) 2 , (17) 

whose proof involves the Gauss theorem on the two surfaces r\ 2 = s surrounding the two singularities. Thus, for 
singular functions, the integral of a gradient is not zero. This indicates that the "ordinary" derivative is not adequate 
for the purpose of application of a "fluid" formalism to point-particles (because the integral of a gradient is always 
zero for fluids). Instead one must generalize the derivative to take into account the singularities, in a way similar to 
the distributional derivative of distribution theory. [[30| 

Another nice connection between the two definitions (3-4) and (3-5) is 

lim Pf / d 3 x SJx - y x )F{x) = (F)x , (18) 

e^O J 

where <5 e denotes the Riesz jl5] delta-function defined for any e > by 

5 e (x) = e<yl ~ € \ x\ e ~ 3 , so that Alaf- 1 = -4ir6 e . (19) 

47T 

Clearly, in the limit e — *■ 0, the Riesz delta-function yields a generalization of the Dirac measure applicable to the 
Hadamard partie finie of a singular function in T . 

To conclude, let us give the result of the application of the general formalism to a binary system of point-particles, 
concerning the total flux (or "gravitational luminosity" C) emitted by the binary in the form of gravitational waves. 
This quantity plays a crucial role in the computation of the orbital phase of the binary, as it evolves with time taking 



4 In addition to (3-3) we assume that the functions F £ T decrease sufficiently rapidly when |x| —* 00, so that all integrals we 
consider are convergent at infinity. 
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into account the loss of energy by gravitational waves. For circular orbits the orbital phase <f> = J udt is obtained 
from the energy balance equation as 

dE n : f udE 

n = ~ c * *=-] — • (20) 

where E denotes the orbital binding energy of the binary in the center of mass frame. Let tn\ and m 2 be the two 
masses, and denote m = vti\ + rri2 and the mass ratio v — \ijm — mim2/m 2 . As a small post-Newtonian parameter 
we define x = (Gmw/c 3 ) 2 / 3 where u> = 2n/P is the orbital frequency and P the period [a; is of order 0(e 2 ) in the 
notation (1-1)]. The present post-Newtonian accuracy for C is 3.5PN except that the contributions proportional to 
the mass ratio v in the 3PN term are not yet under control. These contributions are indicated by 0(y) in the formula 
below; their computation is a work in progress (collaboration with Faye, Iyer and Joguet). We obtain 

r 32c 5 2 sfi ( 124 7 35 \ t 3/2 
r = v x 5 < 1 + — - — - —v x + &rx 3 ' 2 



3G IV 336 12 

44711 9271 65 



v H v 2 ) x 2 



9072 504 18 
8191 535 



-v 



,5/2 



672 24 

/ 6643739519 1712 „ 856, . . 16 



C - ^ ln(16a;) + ^-ir 2 + 0{v) x 3 



V 69854400 105 105 
16285 176419 19897 



504 1512 378 



2 ^^7/2 



V 7TX 



0(x 4 ) . (21) 



The Newtonian approximation was given by Landau and Lifschitz |35| and the classic work of Peters and Mathews |3| 
(recall that the Newtonian term is sufficient to account for the radiation emitted by the binary pulsar PSR1913+16) 



The 1PN term was computed long ago by Wagoner and Will |37| applying the formalism of Epstein and Wagoner, |2C 
and later confirmed |Q applying the present formalism. The 1.5PN term (with Att as a factor) is due to gravitational- 
wave tails , w hich was first obtained by Poisson |3^] using a black-hole perturbation, and it was re-computed by 
Wiseman M] and Blanchet and Schafer |0] within the present formalism. The 2PN term was computed independently 
by Blanchet, Damour and Iyer p2] applying the same formalism, and by Will and Wiseman [|19| applyin g th eir new 



formalism (see Ref. []43j for a summary). The 2.5PN and 3.5PN orders are due to higher-order tail effects. 44 17j] The 
3PN order is very interesting. It involves a logarithmic term (first obtained in this context by Tagoshi and Nakamura, 
0) as well as the Euler constant C — 0.577 ■ • • and a tt 2 . These contributions in the 3PN order, computed in Ref. [fl7|| , 
are due to the so-called tails of tails, that is to say, tails generated by the tails themselves — a purely cubic effect. 
The corrections 0(y) in the 3PN term depend on the full details concerning the regularization outlined previously. 

Using a; as a post-Newtonian parameter is convenient because of its invariant meaning (the same in all coordinate 
systems). Thus the formula (3T0) can be directly compared with the one obtained using a black-hole perturbation 
(which typically describes the background in Schwarzschild coordinates) triggered by the motion of a test particle of 
mass vm in the field of a black-hole of mass m. Using this method Sasaki, Tagoshi and Tanaka pl[-p3[ obtained £ 
to 5.5PN order in the limit v — > 0. In the same limit the formula (3T0) agrees completely with these results to the 
corresponding order. Notably the rational fraction 6643739519/69854400 in the 3PN coefficient comes out exactly 
the same in the post-Newtonian formalism as in the black-hole perturbation approach (collaboration with Iyer and 
Joguet). 
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